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Abstract 

We develop the eikonal perturbation theory to describe the strong-field ionization by finite laser 
pulses. This approach in the first order with respect to the binding potential (the so-called gener¬ 
alized eikonal approximation) avoids a singularity at the potential center. Thus, in contrast to the 
ordinary eikonal approximation, it allows to treat rescattering phenomena in terms of quantum 
trajectories. We demonstrate how the first Born approximation and its domain of validity follow 
from eikonal perturbation theory. Using this approach, we study the coherent diffraction patterns 
in photoelectron energy spectra and their modifications induced by the interaction of photoelec¬ 
trons with the atomic potential. Along with these first results, we discuss the prospects of using 
the generalized eikonal approximation to study strong-field ionization from multi-centered atomic 
systems and to study other strong-field phenomena. 
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I. INTRODUCTION 


Historically, the eikonal approximation was first introduced in relation to light scatter¬ 
ing 0,0. Since then, such theory has been extensively applied not only in optics but also 
in atomic and molecular physics, quantum field theory, high energy physics, etc. The first 
application to modern theories was the analysis of scattering processes of energetic particles, 
as originally proposed by Moliere 0. The idea behind the Moliere approximation was the 
interpretation of particle trajectories as classical, straight line paths, which can be explained 
provided that the de Broglie wavelength of the scattered particle is small compared with 
the size of the scatterer 0-0. The pioneering work of Glauber [73J extended the Moliere 
eikonal to the interaction of fast particles with complex atomic and nuclear systems, and 
established the basis of an exceptionally valuable technique to treat scattering processes at 
high energies. 

After the theoretical basis were set, the eikonal approximation (EA) was successfully 
applied to more complex systems. For example, it was applied for the treatment of the inverse 
bremsstralung heating, which involves the electron scattering by a potential in the presence 
of intense electromagnetic fields 0,0. In the original work of Choudhury and Bakar 0, and 
of Zon Q, the scattering amplitude was obtained by using the EA for the nonrelativistic 
Schrodinger equation and the laser field was considered as a monochromatic plane wave. The 


dipole approximation was extensively used in such calculations. Kristie and Mittleman 10 


have shown that, for ionization in strong laser fields, the electron motion should be treated 
relativistically and the dipole approximation is not good enough to determine correctly the 
electron energy distribution and the transition rate as a function of the laser intensity. In 
order to account for those facts, the systematic eikonal perturbation theory was formulated 
by Kaminski |ll| for the nonrelativistic Schrodinger equation and the relativistic Klein- 
Gordon equation. For nonrelativistic case the dipole approximation for the laser pulse 
was used and applied to free-free transitions, whereas for the relativistic one the finite, in 
general, laser pulses in the plane-wave fronted approximation were considered. Note that 
the approach developed by Kaminski in (T0 was based on the proper-time method, which 
for propagators in quantum theories was proposed by Fock 00, and further developed by 
Schwinger [lO, 10 , and by Fradkin and co-workers 05 - 17 ] (see also the review flO]). 

It has been demonstrated that, for potential scattering in the presence of strong fields, the 
EA has certain limitations. Specifically, such a theory cannot be applied in regions far from 
the potential interaction zone [10. Other considerations restrict the EA just to moderate 


laser intensities 20]. The so-called generalized eikonal approximation (GEA), which is the 


first order term of the eikonal perturbation theory 0 . overcame part of the mentioned 
problems by including certain quantum properties of the system in the eikonal limit. The 
inclusion of such terms extended the range of applicability of the theory to large distances 
from the interaction site 19] allowing a proper treatment of the scattering amplitude. The 
GEA was further extended to solving the relativistic Dirac equation in order to account 


20 


approximately for the electron spin effects 

In Refs. [2ll . 22], the so-called eikonal-Volkov approximation was introduced to treat 


strong-field ionization of atoms. This consisted in including the laser field in full extent by 
means of the Volkov wave function, together with the ionic potential treated within the EA 

23H29J). This approach was further developed by Smirnova et 
It included the dynamical analysis of the elec¬ 


tee, for instance, also Refs. 
al. 00, 0 to describe molecular ionization. 


trons according to complex trajectories, which is a powerful tool for the overall understanding 
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of the process (for complex trajectory method, see, for instance, recent reviews [32, 33j). 

Another way to treat the photoionization by intense laser fields is the so-called strong field 
app roximation (SFA), sometimes recognized as the Keldysh-Faisal-Reiss or KFR theory [ 34 ] 


361 ]. The principal idea of the SFA approach is that, after the photoelectron appears in the 
continuum, it is treated as a free particle interacting with the laser held, ignoring any 
interaction with the nuclear potential during its excursion. Even though the SFA theory is 
known as one of the most fruitful analytical approaches in strong-held physics, its application 
can be justihed for short range potentials only, not for the Coulomb type potentials present 
during ionization of neutral atoms or positively charged ions (see, also the most recent review 
by Popruzhenko j33|). The SFA theory gives, in general, a very good qualitative agreement 
with experimental and numerical results for laser light of moderate intensities in interaction 
with a short range potential 0] but it fails to account for several well-documented spectral 
characteristics, including the wrong predictions obtained for the total ionization rate in the 
static held limit |38 . 

Various methods have been proposed for taking the Coulomb interaction into account 
within the SFA. Except of the aforementioned eikonal-Volkov approach 2lH3ll. |39fl, it has 


been done, for instance, by means of the Coulomb-Volkov anzatz j40|-45j which accounts 
for the asymptotic phase of the atomic field-free wave function. In this context also, the 
saddle point method was reformulated in terms of quantum trajectories by Popov and co¬ 
workers 46-49], by Gribakin and Kuchiev 0, and reviewed recently by Popruzhenko (0 . 
The quantum trajectories together with a detailed analysis of the saddle point equations, 
have shown to make considerable improvements in the SFA to include the Coulomb inter¬ 
actions [51j. In our paper, we use the quantum trajectories as presented in Ref. J38J. Note 
that the extension of this method to the Dirac equation has been also reported [521 ]. 

In the present paper we propose a generalization of the eikonal approximation, in order 
to analyze the photoionization of atoms or ions by short laser pulses. This approach avoids a 
singularity at the center of the binding potential and, in contrast to the EA, can be analyzed 
in terms of quantum trajectories. In Sec. [TT1 we set the physical and mathematical basis 
of the GEA and derive an analytical expression for the ionization amplitude including the 
binding potential interaction. We show how the EA is obtained as a special case of the 
GEA for short times or for large distances from the potential center. In the same way, 
we demonstrate that the first Born approximation can be directly derived from our more 
general approach. Sec. IHIl is devoted to the numerical analysis of the photoelectron spectra 
generated by hydrogen interacting with an intense short pulse. Using both the Keldysh 
theory and the GEA theory, we discuss the coherent diffraction pattern in the photoelectron 
energy spectra and the role played by the atomic potential. We present our concluding 
remarks in Sec. m 

II. THEORETICAL FORMULATION 

For a time-dependent problem described by the Hamiltonian H(t), which for our further 
purpose we separate into two parts, 

H(t) = H 1 (t) + H 2 (t), (1) 

the time-evolution operator U(f, t') satisfies the following Schrodinger equation, 




( 2 ) 
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with the initial condition, = I. Here, I is the identity operator. The solution to the 

above equation, which incorporates also the initial condition, can be written in the form 

U(f, t') = Texp^—i J drH(r)^, (3) 

where T is the time-ordering operator [b, 53}. Let us also introduce operators Ui (£,£') and 
U 2 (£, £') which determine the time-evolution governed by the Hamiltonians H^f) and H 2 (£), 
respectively. In other words, for each of them it happens that 

i^lb(M') = £'), U *(£',£') = i, (4) 

and 

U i(t,t') = texp^-i j drHj(r)^, (5) 

for i = 1,2. It can be shown that the propagator U(£, £') given by Eq. (I3j), fulfills the integral 
Lippmann-Schwinger equation such that 

U(M') = Ui (£,£') - i [ dr U(£, t)H 2 (t)Ui(t, t'). (6) 

Jt' 

Since now on, we will assume that the Hamiltonian H^i) is independent of time, H(f) = H l5 
whereas the Hamiltonian H 2 (t) varies with time in the interval when t G [0,T] and is zero 
otherwise, 

H 2 (f) =0 for t < 0 and t > T. (7) 

Our aim is to calculate the probability amplitude for a system governed by the Hamiltonian 
H(f) to make a transition from the initial state |^) to the final state |^f)- We assume that 
these are stationary states of the Hamiltonian Hi, 

HijV’f.i) = E , f , i |'0f,i), (8) 

which are orthogonal = 0. For completeness, we note that their time-evolution is 

given by 

|Vk,i (t)) = Ui(t, O)|-0f,i) = e _lE W|V'gi}. (9) 

The aforementioned probability amplitude calculated at time T is 

A,i(T) = (^f(T)|U(T,0)|^(0)). (10) 

It follows from Eq. ((HD that this quantity can be rewritten as 

Ai (T) = (MTM(T))-i f T dt 1 ('0 f (T)|U(T, f')H 2 (t')|-0i (;£')), (11) 

J o 

where the first term vanishes due to the orthogonality of the initial and final states. 
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A. Ionization probability amplitude 


The above theory can be conveniently applied to describe short-pulse ionization processes. 
For this purpose, let us specify that is the atomic Hamiltonian, 

H ^^ 2 + v ^ (12) 

whereas H 2 describes the coupling with the laser field which, in the length gauge, has the 
form 

H 2 = —eS(t) ■ r. (13) 

At this point, we also specify that the vector potential, which describes a finite laser pulse, 
A(t) depends on time for t e [0, T] and is 0 otherwise. It is related to the electric field £{t) 
such that A(t) = — dr£ (r). The initial state is the atomic ground state of energy E 0 , 
-00 (n), which evolves in time according to 

(r|0i(t)) = e“ lE °Vo(r). (14) 

The final state is the scattering state 0p 1 ( r ) which describes a particle of momentum p, 

2 

(r|0f(t)) = exp^i^-^0H(r). (15) 

Next, let us define the retarded propagator in the length gauge K L (r,t ; r',t') which relates 
to the quantum-mechanical dynamics governed by the Hamiltonian H(t), 

K L (r,t-r',t') = 9(t - t')(r\U(t,t')\r'), (16) 

with the initial condition that 


K L {r,t + 0; r',t) = 8{r - r'). (17) 

With these definitions, the transition probability amplitude (ITTli for ionization becomes 

/ t 2 r r 

dt'exp Sr / Sr'il>^(r)K L {r,T\r',t')(-e£(t')-r')ip 0 (r'). 

(18) 

Therefore, for our further analysis it is necessary to derive an explicit form of the propagator 
Ki(r,t-, for an electron under a simultaneous action of the laser field and the external 

potential. This will be done in the next Section, using eikonal perturbation theory El). 


B. Retarded propagator under the GEA 

While our ultimate goal is to derive the retarded propagator in the length gauge, we start 
with the velocity gauge and the respective propagator Ky(r , t; r', t'). For the time being we 
consider the most general case in which the scalar potential also depends on time, V(r,t). 
In the present case, Ky(r, f; r ', t') satisfies Eq. (TT6]l but with the Hamiltonian written in the 
velocity gauge, 

(- 1 J 7 ~ - eA ^')\ 2 - V(r',t')^K v {rit]r',t') = i S(t - t')5(r - r'). (19) 


5 


Then, one can look for the propagator Ky{r,t ; r',t') using the Fock-Schwinger proper-time 


representation 11 


. , A r °°, fdnd 3 k 

K v (r,t,r,t)=J o d s J-^e W 


k 2 

-i Q(t — t') + ik ■ (r — r') + is(fi — —-f ic) 


+is) + i Xk(r',t',s) 


( 20 ) 


with s being the proper time. Here, we have introduced unknown functions s) and 

Xk( r 'i t 1 , s). These functions must satisfy the following conditions 


, s) = 0 when A = 0, 
Xk{r', t', s) = 0 when V = 0, 


( 21 ) 

( 22 ) 


in which case the propagator (120]) becomes a free particle propagator 11|, [18]. When sub¬ 
stituting the formula (120]) into Eq. dl9]) we arrive at 


dfld 3 k 


(2tt) 


-i£l(t—t')+ik-(r—r') 


ds 


i 2 2 

n -- + —fc • A(t ') - — A 2 (t') - V(r', t') (23) 

2m m K J 2m y J v ’ ' K J 


+d t '$k + d t 'Xk + —{k - eA(t')) • V'y fc + ^—A'x k 
m 2m 

k 2 


2m 


(v'W 


x exp 


is ( fl - --b is ) + i$fc + i Xk 


= i 5(t — t')8(r — r'). 


If $fe(t',0) = 0 and x fc (r',f', 0) = 0, the equation above is fulfilled provided that 


ds d s exp 


k 2 

is (H - — + E ) + i$fe + i Xk 


= i , 


(24) 


As a result, we obtain a partial differential equation for functions <&*.(£', s) and Xfc(r',f',s), 
which then can be separated into two independent equations: 




(25) 


^p~§~) Xk ( r> ' ^ s ) = “ ~ [ fc “ eA (0] ' V'x fe (r', t', s) + H(r', t') (26) 


+ 7 ^(V , y fc (r',t',s)) 2 - ^-A'y fc (r',t',s). 


When solving these equations one has to remember of the above initial conditions for s = 0 
and about Eqs. (l2Tj) and (122]) . It is rather straightforward to solve Eq. (125]) . On the other 
hand, Eq. (126]) is a nonlinear second order differential equation for Xk( r ',t'i s), which can 
be solved explicitly only for particular potentials, for instance, for the harmonic oscillator. 
Since now on we will call Xki'f' 1 , t 1 , s) the eikonal. 

It follows from Eq. (1251) that 


$k(t',s)= ( dr — A(t) • [k — ^A(t)\. 

J t i m L 2 J 


(27) 
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Taking this into account we conclude that in the absence of external potential (in which 
case Xk = 0), the exact propagator (120|) reduces to the Volkov propagator. In order to solve 
Eq. (|26|) perturbatively, let us first rewrite this equation such that 

(J 7 “ Ts + ~ • v ' + i A ') Xfe(r '’ ( '' s) = Wk(r ’’ t '' s) ’ (28) 

where we have defined, 

W k (r’,t',s) = + -L(V' Xfc (r',f', S )) 2 . (29) 

Note that in the case when V(r\ t') = 0, it follows from Eqs. (l22l) and (|29l) that W k (r', t', s ) = 
0. Now, with the help of the Fourier transforms, 

»(*•', t',s) = (30) 

W k (r',t', a)= I s), (31) 

we replace Eq. (12%)) by 


d d 

-7W + 


. dt' ds m 
We look for the solution of this equation in the form, 


2 

■ [k - eA(t')} - ^-jx fc (/^,f',s) = W k (K,t',s). 


(32) 


Xk{v,t',s) = exp i J dr (£ • [k-eA(r)\ - —) x! k {n,t',s). (33) 

By putting this solution into Eq. (]32l) . we find out that a new function x' k (n,t\ s ) satisfies 
the following equation, 


rt'+s 


- dt' ds 
which leads to 




^--^)x' k (^ t ’,s) = W k (ii,t',s)ex.p\-i I d r(£-[fc-eA(r)] 


rt'+s _ r rt'+s 

dcr W k (/J,, a, t' + s — a) exp —i / dr( — [k — eA[r)\ — 


m 


2m J. 


(34) 


(35) 


From here, it also follows that 


Xk(p,t',s) = - exp 




daW k (^,a, t' + s 


x exp 


rt'+s 


■/ da '(+,- [k - eAia "> ] --L)y 


Before we proceed further, we introduce new quantities 

1 

a k (t)=— d r[fe-eA(r)], 

Tfl ./n 


R k {r' : t', a) =r' + a k (a) - a k (t ') = r' + — f d r[k - eA(r)]. 

m Jv 


(36) 


(37) 

(38) 
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Then, Eq. (I36|) can be rewritten in a more compact form, 

rt'+s 

Xk(vJ,s) = - \ 


daW*,^, a,f' + s-a)exp(i/i- [a k (a) - a k {t’)\ - i|^-(a - t ')), (39) 


which, after substituting into Eq. (l30]h gives 

rt'+s r 


Xk{r',t',s) = 


x 


j da J d 3 p W k (p, a, t' + s — a) 
^jexp(i^ ■ [i4(r',f» - p] - ‘£^(a - t!) 


(40) 


Here, we recognize that the integral over p is the Fresnel integral, and so it can be performed 
exactly. In doing so, we arrive at the following expression for the eikonal, 


»(r', t\ .)=-/ da j dV ( M( ”_ t0 )' /2 [R k (r', t\ a) - pf) 

ft '+s 


xW k (p,a,t' + s — a) = 


da V eS (R k (r', t', a), t', a), 


(41) 


which implicitly defines an effective potential V e s. This is the starting point for the eikonal 
perturbation theory and for the GEA. 

For completeness, let us go back to Eq. (120|) and perform the respective integrals over fl 
and s. As a result, we obtain the integral representation of the retarded propagator in the 
velocity gauge such that 


K v (r,f,r',t ') = 


d 3 k 

(2^)3 


exp 


i k ■ (r — r') — i(t — t ') 


2m 


(42) 


+i t-t') + i Xk(r', tt - t') 
with the functions & k (t', s ) and Xfe(r', t', s ) defined by Eqs. ()27li and (ITTTh respectively. Since 

-dR k (r',t',a)' 2 


2 m y J ; 2 


da 


da 


we can rewrite Eq. (ITZji such that 

d 3 k 


K v (r, t;r ,t) = 


(2jt) 


—i / da 


exp|ifc • (r — r') 

m fdR k (r ', f', a) 
da 


(43) 


(44) 


+ V eS (R k (r' a),t', a) 


On the other hand, one can show that in the length gauge the retarded propagator has the 
following form, 


K L (r , t] r',t') = 


d 3 fc 

(2vr) 3 


exp 


i(fc — eA(t )) • r — i(k — eA(t ')) ■ r' — i (t — t ')—— 

2m 


+id>fc(t', t-t') + iy fc (r', t', t - t’) 


(45) 
















with exactly same functions & k (t',s) and Xk(r',f , s) as before. This, in turn, leads to 


, , f d 3 k f. 

K L {r,t\r ,t) = J ^-^expji m 


dR k (r',t',t) . dR k (r',t',t ) , 


dt 


■ r + \m 


dt' 


■ r 


(46) 


—i / da 


m /dR k (r’, t ', a) 


L 2 


da 


+ V eS (R k (r',t\a),t\a) 


and is the most general form of the propagator defining the ionization amplitude in Eq. (TT7T) . 

In the first order approximation with respect to the potential, W k (r', t', s ) can be approx¬ 
imated by V(r',t') which is the essence of the GEA. Hence, we can write that, in the first 
approximation, the eikonal (1411) becomes 


XfcV'^s) = - / daV^ } (R k (r , ,lf,a),t , ,a), 


ft' +S 


l 1 ) t (nr,' +' +' 


(47) 


where 


Ve&{ R k( r ',t',a),t',a) = I d 3 p( 


m \ 3 / 2 


27ri(cr — t’) 


(48) 


We will also use the respective notation for the propagators Ky * (r, t; r', t') and K ( r ^ (r, t ; r', t') 
Note that for the static potential, when V(p, a) = V(p), depends explicitly only on 
R k and a — t'. 

At this point, let us investigate the limit of short time intervals, i.e., when a ~ t': in this 
context the short-time limit is equivalent with the classical one, when h —>■ 0. To this end, 
in Eq. (148|) we make use of the following model of the delta function 


5(r) = lim 

e ->0 


exp 



(27rie) 3 / 2 ’ 


(49) 


which leads to 5(R k (r',t',a ) — p). Performing the remaining spatial integral in (148|) . we 
arrive at the conclusion that 


Vrt( R k(r',t',cr),t',a) ss V(R k (r',t',a) t a), 


(50) 


which holds in the limit when a ~ t'. This result is in full agreement with the Dirac 


conjecture [54j, |55j that in the short-time limit (and only in this limit, i.e., when the quantum 
spreading of the electron wave packet is negligible) the propagator is proportional to e l5cl , 
where S c \ is the classical action. This afterward has lead Feynman to the path integrals jstij . 

In closing this Section, let us also note that the integral equation for the eikonal (fXT|) 
allows to construct a series expansion with respect to the potential - the eikonal perturbation 
theory. Note that such theory and the GEA, as its first order term, was proposed in Ref. 13 
for both the nonrelativistic Schrodinger and relativistic Klein-Gordon equations, and it 
was applied to free-free transitions. In this paper we further extend it for nonrelativistic 
ionization processes. 
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C. Generalized eikonal for the Coulomb potential 


For the Coulomb potential describing the interaction of an electron of charge e < 0 and a 
nucleus of charge — Ze, where Z — 1,2,... is the atomic number, V(r,t ) = V(r) = — Zac/r. 
Here, a = e 2 /(47T e 0 c) is the hne structure constant; in the atomic units used in our numerical 
analysis ac—1. The eikonal defined by Eq. (PT7T) becomes 




Zac 

x-exp 

P 


1 m 


2 (a- t') 


[■ Rk(rt',a) 



(51) 


where the integral over p can be performed exactly. This leads to 



rt'+s 


Zac 


da ■ 


\R k {r',t',a) 


r erf 


m 


2 i (a — t') 


| R k (r',t',a) 


(52) 


where erf(z) is the error function. The commonly used eikonal (see, for instance, in 0 , 0 , 
[111 . Ill'] and references therein) which, contrary to the above approximation, is singular for 
the Coulomb potential, 




Xfc,original^ ? t ? <§) ZOLC I d(J 


R k (r’, <r)| ’ 


(53) 


is recovered, if 

777 

^—^|i4(r / ,t , ,a)|»l. (54) 

Note that this condition is in agreement with the short-time interval approximation, a ~ 
t', which has been discussed above. However, for a ^ t', the aforementioned condition 
requires that \R k (r', If, a)| 0. In other words, the approach proposed in this paper is 

also applicable to cases when the electron trajectory can return back to the origin of the 
Coulomb potential. For this reason, it does not lead to problems already mentioned in 
Ref. [38j. In our case: (i) the integral in for a close to t' converges (for the integrand 
we meet only the integrable (a — t')~ R 2 singularity), and (ii) the generalized eikonal is not 
singular for the zeroth-order trajectories that may revisit the nucleus in real time. It is 
worth mentioning that the condition (l54j) is fulfilled for large distances from the Coulomb 
center. This means that the original and generalized eikonals coincide with each other 
not only for short times, but also at distant points in space. Hence, if the electron wave 
packets or quantum trajectories are far away from the center during the time evolution both 
approximations should give similar results. This is usually the case if the final kinetic energy 
of photoelectrons is much larger than 3U p , where U p is the ponderomotive energy defined 
below [Eq. (1771) ]. This will be demonstrated later on. 

In closing this Section we note that similar close expressions for the eikonal X k \ r \t\s) 
can be also derived for other potentials such as the Yukawa, Gaussian or multi-center 
Coulomb potentials. These cases will be studied in due course. 
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D. Ionization probability amplitude in the GEA 


In the first order eikonal approximation, the ionization probability amplitude in the length 
gauge equals 

rT 2 


2 

A {1 \p) = - i / df'exp(i^-T - i E 0 t') 

Jo 2m 

x [ d 3 r [ d 3 r'^~ ) *(r)K { i\r,T-,r',t')(-e£ (0 ■ OV’oIOj 


(55) 


with 




d 3 /c 

(2vr) 3 


exp 


i k ■ r — i(fc — eA(t')) ■ r' — i(T — f') 


2m 




(56) 


Having this in mind and performing the integral over r in Eq. (loo]) , we arrive at 
rT - 2 r r d 3 k 


= dt / exp(i^—T — iE 0 t') j d 3 r' J ) *{k)(-e£(t') ■ r')^ 0 (r') (57) 


x exp 


-i (k - eA{t')) ■ r' - i (T + i§ k {t',T - t') + ix k \r',t', T - t') 


In the following, we assume that the Fourier transform of the final scattering state, 'tpp \k), 
is centered around k ~ p. Essentially, this corresponds to a plane wave approximation in 
the final electron state which is the common approximation applied in strong-field physics. 
In other words |Jjp ^*(0)|, which emerges under this approximation in Eq. (1571) . contributes 
only to the normalization of the plane wave. Therefore, we can disregard this multiplication 
factor remembering that the density of final electron states per unit volume equals d 3 p/ (27t) 3 . 
This way, the probability amplitude of ionization (157]) becomes 


A {1} (p) = - i / dt ' exp i(^- E 0 )t' + i$ p (t',T - t') 

In L Zm 


(58) 


x J d 3 r'(—e£(t') ■ r^ipoir') exp — i(p - eA(t')) ■ r' + iXp\r\ t', T - t') , 
where <hp(f',T — t') and Xp\ r 'A'->T — t') are defined by Eqs. (]27j) and (]52]h respectively. 


E. Limit of the Born approximation 

In order to establish relations between the GEA and the Born approximation, we go back 
to Eq. dl8]) . As in the previous Section, we assume that the final scattering state is described 
by the plane wave of momentum p, i.e., cs e ipr . The retarded Volkov propagator in 

the length gauge, r*', t'), can be obtained from Eq. (145]) by neglecting Xk • Hence, 

d 3 k 




(2 Jr) 


exp 


i (k — eA(t )) • r 


(59) 


rt i 

i(fc — eA(t')) ■ r' — i / dr—(fc-eA(r)) 2 
7 / 2 m . 
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The wave function ipp\r', f'), defined by the integral 
^(o)*(r',f')= [ d 3 re- ip r K^\r,T-r',t') 


( 60 ) 


= exp 


— i(p -» eA(t')) ■ r' — i / dr-—(p — eA(r))‘ 


is the Volkov solution of the Schrodinger equation in the length gauge. One can check that 
it fulfills the boundary condition 

*Pp\ r ' ,T) — e‘ pr ‘, (61) 

as for t' ^ T the action of the laser pulse vanishes. 

The exact retarded propagator, K L (r, t; r', £'), satisfies the Lippmann-Schwinger equation 
in accordance with Eq. ([6]), 


K L {r , t; r\ t’) =K^ ] (r, t; r', t') 


(62) 


-i J d r J d i yKf ] (r,t]y,T)V(y,T)K L (y,T\r',t'), 

which allows to split the exact probability amplitude [cf. Eq. (U5|) ] into two terms, 

A(p) = A k {p) + A resc (p)> 

where 

A K (pi = -i / dp exp /dVVi<” , *(r',P)(-ef (if) • r')p> 0 (r') 


(63) 


(64) 


is the Keldysh amplitude, and 
Aesc(p) = ~ dt' exp 


T 2 

1 - 1 — 1 Lnt 

L 2m 


dV / dr / dV0( o) *(t/,T)V(t/,T) (65) 


xK L (y,T]r',t')(-e£(t') ■ O^o(r') 


is the exact rescattering amplitude. In the first Born approximation [by the Born approxima¬ 
tion we understand the expansion with respect to the potential V(r, t)] the exact propagator 
in the above equation is replaced by the Volkov propagator in order to get the amplitude that 
describes the rescattering of electrons after ionization. After some algebraic manipulations 
one can show that 


Kf[y,r t r\t’ 


m 


i>f'(y,r) K27ri (t-0 


3/2 


( 66 ) 


x exp 


i nn 


L2(r - 1 ') 


Itpir,If. T) -y 


Inserting this identity into (1651) . we can write down that the rescattering amplitude in the 
first Born approximation is A reS c(p) ~ A&\(p), where 


Abi(p) = f dV f dt' exp i^—T — iE 0 t' 


(67) 
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and Xp\ r 'i s ) is defined by Eqs. (I47|) and (148]) . It follows from Eq. (158|) that the total 
probability amplitude in the first Born approximation, Ak(p) +*-4bi(p), is exactly recovered 
from our first order GEA provided that | Xp\ r ',t',T — t')\ <C 1. This means that one can 
use the Taylor expansion 






r'i t', T — t') 


( 68 ) 


Such an agreement is not achievable within the EA. This shows that the first Born ap¬ 
proximation, which has been extensively used in the analysis of rescattering processes in 
ionization [32|, is the limiting case of the generalized eikonal expansion. 


III. COMBS IN THE PHOTOELECTRON SPECTRUM 
A. Laser pulse and its characteristics 

In the above formulation we have assumed that a finite laser pulse lasts for time T and, 
therefore, it is described by the electric held £{t) which vanishes for t < 0 and t > T. The 
pulse duration T defines the fundamental frequency of held oscillations, c o = 2n/T. One can 
also introduce the held phase, 0 = u>t, which allows to rewrite the above condition such that 
£(0) vanishes for 0 < 0 and 0 > 27r. We assume that the driving pulse is linearly polarized 
along the z-axis. In the dipole approximation, the laser held is described by the electric held 
vector 

£{4>)=£o fs(<l>)ez, ( 69 ) 

where £q is related to the amplitude of held oscillations. Here, the shape function is 

adjusted such that fg (0) = 0 for 0 < 0 and 0 > 27 t, and it has to satisfy the condition [32] 

/ /e(0) d0 = O. (70) 

Jo 

This condition is fulfilled provided that the shape function has the following Fourier decom¬ 
position, 

N 0 

feW = Y! ( 71 ) 

N=-N 0 

where Yl' means that the zeroth Fourier component is excluded from the sum, JV^O. Since 
/f(0) is a real function, we also require that = £-n- This expansion allows us to define 
the average intensity carried out by the laser pulse, 

1 f 27T 

I = ice o£ 2 ) = — J ce o £ 2 (0)d0. (72) 

Namely, 

N 0 

I = 2ce 0 £^\£ N \*. (73) 

JV=1 

Note that this definition, even though derived for a finite laser pulse, is consistent with the 
monochromatic plane wave approximation. In the latter case, taking the electric field of the 
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form £(4>) = £osin (<f>)e z , we have £n = =Fi/2 for N = ±1. Hence, it follows from Eq. (175]) 
that the averaged intensity of the monochromatic plane wave equals / = ceqEq 2 /2. 

Let us now define the shape function for the vector potential A(0), 


/u(</>) = - / fe(<p)d<P, 


(74) 


which leads to 


AW = —fA{4>)e z . 

CO 


(75) 


Using Eq. (1711) one can derive the corresponding Fourier decomposition of the shape function 
/u(0) [Eq. C3D] 


N 0 

/a(0) = A) + A N e 


-iN<j> 


(76) 


N=-N 0 


N 0 


where A 0 = —2 ReA^ and A N = —i£ N /N assuming that N ^ 0. Since we have imposed 


N =1 


the condition A(0) = A(27t) = 0, the vector potential has a constant and an oscillatory 
contributions, AW = A const + A osc (</>). It is the oscillatory contribution to the vector 
potential, A osc (0), which describes the quiver motion of free electrons in the laser field. The 
ponderomotive energy of such a motion can be defined as 


e 2 A 2 

U p = ( osc 


r- 2 tt 


2m ^ 27r 




(77) 


Applying here the series expansion (1761) . we End out that 


e 2c2 N ° I c 12 

e c o Rjv| 

moo 2 jV 2 

N=1 


(78) 


Again, for the monochromatic plane wave, we obtain from Eq. (1751) that U p = e 2 Eq /(4mo; 2 ). 
This is the well-known formula for the ponderomotive energy of a free electron driven by 
the monochromatic plane wave. 

For our further purpose, we introduce also the vector function 


where 


with 


«(</>) 




faW e z, 


r4> 

faiW) / f a{W) d^2 fa,0 T fa, 10 T fa, osc(0); 

JO 


(79) 

(80) 


fa, 1 : 
fa, osc(0) 


(81) 

(82) 

(83) 
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This function will be used in Sec. IIII C II 

In this paper, we consider the laser field described by the shape function 




sin 2 (iV rep f) sin(iV rep 0), 

0, 


0 e [o, 2 tt], 
otherwise. 


(84) 


Such a laser field consists of N rep single-cycle pulses with no time delay in-between (where 
iV re p = 1, 2, 3,...). Introducing the laser frequency cul = N rep uj, we can represent Eq. ()84l) as 


fe(t) 


\ sin(cu L t) 

0, 


| sin(2w L t), te[0,T], 
otherwise. 


(85) 


This clearly shows that the laser field (I87|) can be composed out of two harmonics, i.e., (Ul 
and 2 cul. For our choice of the shape function, the only nonzero coefficients in its Fourier 
expansion (171]) are £± 7 v rep = Ti/4 and £± 2 iv re p = ±i/8. Therefore, according to Eq. (173|) . the 
averaged intensity carried out by the laser pulse (184|) is 

1 = ^2 C£ °^' ( 86 ) 

Moreover, it follows from Eq. d78|) that the ponderomotive energy associated with the quiver 
motion of an electron in such a held equals 


rr 17 e 2 £g 

p 256 mul 


(87) 


Further we shall assume that cjl equals the frequency of the Ti-Sapphire laser, cul = 1.55eV, 
while the averaged intensity of the pulse is / = 3.125 x 10 13 W/cm 2 . Thus, the ponderomotive 
energy of the electron oscillating in the laser pulse (1841) equals U p = 1.024 o ; l - 

Below we shall analyze the energy spectra of photoelectrons ionized from a hydrogen-like 
atom by the laser held (IMP with different N rep . This means that while changing the number 
of pulse repetitions, A^ rep , the time duration of the entire sequence of pulses, T = 27rfV rep /a;L, 
will change as well. First, we will present the respective results based on the Keldysh theory. 


B. Combs in the Keldysh theory 


In the Keldysh approximation, the amplitude of ionization Ak{p) is given by Eq. 
This formula can be written explicitly in the form, 


A K (p) = -iexp{i(f^ 


EAT- i 


dr 


'^ p - eAi E 2 - E '‘ 


-i(p — eA(i')) ■ r' + i j t^(p - hA(t)) 2 - E 0 |(-e£((') • (88) 


x exp< —1 


Introducing here the phase of the laser pulse, 0, and defining the following quantities: 


q(0) =p-eA{<j>), 


G{p, 0) = - / d0' 

OJ .In 


2 m 


En 


d t' J d 3 r' 


r')ipo(r'). 

(88) 

quantities: 



(89) 


(90) 
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we can rewrite Ak{p) as 

pi'Mp) r 2n r 

A K {p) = —i——— J d (j>e iG{p ^ J dVe- iqW ' r (-e£(0) • r')^ Q {r'). (91) 

Here, we have also introduced the following abbreviation, $ 0 (p) = — E 0 )T — G(p,2n). 

For a hydrogen-like atom in the ground state we have 


M r ) = M -e 


—A r 


En = — 


2m ’ 


(92) 


where A = (Za 0 ) 1 and a 0 is the Bohr radius. In this case the space integral in Eq. (19T1) can 
be performed exactly. As a result, we obtain 




where, according to Eq. (l90lh 


d 


G'(p.0) = ^-G(p,0) = -[^h) 

0(p wL 2m 

Q 2 \ 

G"(p,4>) = =- ^e£((p) ■ q{<f>). 

o<p z mw 


-En 


(93) 


(94) 

(95) 


Inserting Eq. (1931) into Eq. (T9TT) . we obtain that the probability amplitude of ionization of a 
hydrogen-like atom in the ground state, within the framework of the Keldysh theory, equals 


Ak(p) 


4A 2 


\/A7r 
m 2 u.i 2 


gi^o(p) 



G"CM) 
[G'(P, 0)] 3 


gi G{p,<t>) 


(96) 


To calculate the total probability of ionization, one has to integrate \Ak{p) ' over the density 
of final electron states, d 3 p/(27r) 3 . In doing so, we arrive at the following formula, 


V K = 


d 3 p 

(2^j 3 


\Ak ( P )\ 2 = 


m 


(2vr) 


dH p / dE p \p\\A K (p)\ 2 , 


(97) 


where E p = p 2 / (2m) is the electron kinetic energy and dh2 p is the electron differential solid 
angle. The integral remaining in Eq. (1961) can be performed numerically. Before proceeding 
with numerical calculations, let us derive also an approximate formula for the probability 
amplitude (jU5|) using the saddle point method. This analytic approach will help us later on 
to interpret our numerical results. 


1. Singular saddle point approximation 


In light of Eq. (19H we consider the integral, 


Z v 



[G'(p, 0)] I/ 


gk G(p,4>) 


(98) 
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where f/(0) is a regular function of 0 and £ > 0 is a large parameter. The saddle points 0, 
are the solutions of the equation, 


(99) 

While evaluating the integral in Eq. (1961) at 0 = 0 S , we encounter the problem in the 
denominator which contains G'(p, 0). Therefore, to avoid singularities at the saddle points, 
in Eq. (198|) we substitute, 


[G'(p, - 




dr] T] 


^-1 e ~G'(p,cj>)r) 


( 100 ) 


Then, an ordinary saddle point method can be successfully applied which consists in replac¬ 
ing the function 77(0) by its value at the saddle point 0 S and, also, in replacing the argument 
of the exponent by the first nonvanishing terms arising from the Taylor expansion around 
0s 33), [50|. After performing the remaining integrals we obtain that 


T ^ 

- Lf l/ 


^(-d/2 

< i^+i \ / v 1 


i( , + i W4 i(^i 

\G"{pM v +W 


2b'-i)/2r(£±i) 

where the sum is over such saddle points that satisfy the conditions: 

Im[G(p, 0s)] > 0, Im [G"{p, <j> a )\ > 0. 


( 101 ) 


( 102 ) 


As it will follow shortly, these conditions are compatible with the requirement that lm(0 s ) > 

0. 

Going back to Eq. (1961) and making use of Eq. (11011) for u = 3 we obtain that, under the 
saddle point approximation, the probability amplitude of photoionization from the ground 
state of a hydrogen-like atom equals 


^(saddle) _ 


'■K 


(P) = -2 


A / 7tA \ 2 


,i^o(p) 


7T 


E 


gi G(p,<t>s) 

G"(p, 4>,) 


(103) 


This dehnes the total probability of ionization, in accordance with Eq. (197)1 . Because the 
contributing saddle points have to satisfy the conditions (j 102)1 . their careful analysis is 
necessary. 

For our choice of the pulse shape (J84)), there are in general 8N rep solutions of the equa¬ 
tion 1)99)1 . However, a half of them does not fulfill the conditions (110211 . Among the remaining 
solutions we can distinguish two groups of solutions (2N Tep points each) with the exact same 
positive imaginary parts for 0 S , G(p,(j) s ), and G"(p,(p s ). To illustrate this, we consider a 
Ti:Sapphire laser (cul = 1-55 eV) producing a held composed out of three single-cycle pulses 
(We p = 3), with the electric held described by Eqs. (ES)) and (j83ft . We choose the averaged 
intensity in the pulse / = 3.125 x 10 13 W/cm 2 . Due to cylindrical symmetry of our problem, 
the positions of saddle points do not depend on the azimuthal angle of ionized photoelec¬ 
trons, just on their polar angle 9 p . Here, we choose 9 p = 0.27T. In Fig. Q3 we plot the real 
(upper panel) and imaginary (lower panel) parts of the solutions to Eq. (1991) which obey the 
conditions (11021) . For N rep = 3, we observe 12 such saddle points. The saddle points which 
are represented by the dashed line have larger imaginary parts than the saddle points which 
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FIG. 1. (Color online) Positions of the saddle points <f> s as a function of the kinetic energy of 
electrons E p , calculated from Eq. (1991) . Only these saddle points which satisfy (11021) are plotted. 
The saddle points with the same imaginary part are marked either as solid or dashed lines. The 
parameters of the driving laser field [described by Eqs. (lUH and AMD ] are wl = 1-55 eV, N rep = 3, 
and / = 3.125 x 10 13 W/cm 2 . The final electrons are detected asymptotically at the polar angle 
dp = 0.27r. 


are represented by the solid line. For this reason, the contribution of the former points to 
the sum in (1101 jl is marginally small and can be disregarded in our further analysis. 

In Fig. El we present the shape functions /a(</>), and f a (4>) for N vep = 3. The 

vertical lines represent the real parts of saddle points, Ke<f) s , for the electron kinetic energy 
E p = 12.22cul ~ YZUp. The thin black lines correspond to those saddle points that do not 
contribute much to the probability amplitude of ionization ( 11031 ) . This is not surprising as 
their real parts correspond to the nearly zero value of the electric field. On contrary, the 
remaining vertical lines (thick solid and dashed lines) correspond to the saddle points that 
have to be accounted for in Eq. ( 1103 |) . We see that the important saddle points have their 
real parts which correspond to the nearly extreme values of the electric field. 

In Figs. E] and IH we draw the dependence of the functions G(p, 4> s ) and G"(p, 4> s ) on the 
kinetic energy of photoelectrons E p for those saddle points (j) s that contribute significantly to 
the probability amplitude of ionization (11031) . These saddle points were denoted in Fig. Q] by 
the solid lines. Among these points, we can distinguish between the ones that relate to the 
maxima (solid blue lines) and minima (dashed red lines) of the shape function fe(4>) (see, 
Fig. ED- We denote these saddle points as 4> y Nrep and 4 > Niep , respectively, with £ = 1, 2,..., N rep . 
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FIG. 2. (Color online) Shows the shape function fs(4>) for a triple (lV r ep = 3) laser pulse defined 
by Eq. (151)1 . In the lower panels, the corresponding shape functions f a (</>) and f a (4>) are displayed 
[Eqs. (171)1 and ([50]) . respectively]. We have marked the real parts of saddle points, Re0 s , as vertical 
lines. While the thin black lines correspond to the position of these saddle points which contribute 
very little to the probability amplitude of ionization (11031) . the major contribution there comes 
from the saddle points marked as the thick (both solid and dashed) lines. The positions of Re</> s 
are for E p « 12U P and 6 P = 0.27r. 



R’p/wl 


FIG. 3. (Color online) Shows the dependence of the real (upper panel) and imaginary (lower panel) 
parts of G(p, (f> s ) [Eq. (1901) ] on the photoionized electron kinetic energy, E p , calculated for the same 
parameters as in Fig. |T| Only these saddle points are accounted for which significantly contribute 
to the probability amplitude of ionization (110311 . 
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FIG. 4. (Color online) The same as in Fig. [2] but for the function G"(p, (f> s ). 


They have the same positive imaginary part but their real parts differ such that 


Re< ^Svip = < / > o + ^(^-l), (104) 

1 ’rep 

Re <Si p = -4>o + -^~ L ( 105 ) 

4 V r ep 

Here, 0 O denotes the real part of the first saddle point which gives a significant contribution 
to the probability amplitude. It can be anticipated from the upper panel of Fig. [3] that 

Re[G(p, = G 0 (p) + 2?r (( - l)F(p), (106) 

R*[G(p, 4>$J] = Go(p) + Mt - 1 )F(p). (107) 


Even though it is possible to derive the exact forms of functions Go(p),Go(p), and F(p), 
it is not of a particular interest. As we will show shortly, only the structure of the func¬ 
tions Re[G(p, (pf)^ J] and Re[G(p, )] is important for interpreting the resulting energy 
distributions of photoelectrons. Moreover, it follows from the bottom panel of Fig. [3] that 
Im[G(p,</>55 )] = lm[G(p,<j$ )] = W(p) > 0. Another important observation, based on 
Fig. [4] is that 

Re[G"(p, p )] = —Re[G"(p, (108) 

Im[G"(p, 0 wlep ) 1 = Im[G"(p, $5^)]. (109) 

Having this in mind, we shall denote in the following: Gq(jp) = \G"(p, 0^ jl = |g"(p,$£ji 
and V>C"(P) = arg[G"(p, 0™ <p )l = w ~ ar s[ G "(P' $2.,)]- 

The aforementioned properties of the saddle points along with the discussion of functions 
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FIG. 5. (Color online) Shows the energy spectra of photoelectrons (11111) ionized by the pulse with 
the sin 2 envelope (|84l) . The frequency of the laser field is taken oj l = 1.55 eV and its mean intensity 
is I = 3.125 x 10 13 W/cm 2 . The envelope of the spectra (solid black line) corresponds to a one-cycle 
driving pulse (lV re p = 1). Other curves correspond to a sequence of either two one-cycle (lV re p = 2; 
dashed red line) or three one-cycle (N rep = 3; solid blue line) driving pulses. All results have been 
divided by lV r 2 p . They have also been multiplied by to magnify the main features of 

the distributions. The results in the upper frame have been obtained by performing the integral 
in Eq. (1961) exactly. The results in the lower frame have been obtained using the saddle point 
method (11031) . The major features of the mirror-reflected distributions are the same, except that 
they differ in magnitude. 


G(p, (ft s ) and G"(p , 0 S ) allow us to rewrite Eq. (11031) such that 


A 


(saddle) _ /A /2 ttA\ 2 i[<&[)(p)+7r(jVreD _ 1)F(p)+ l Gn(p)+ lg o(p)1 

K V vr V mu> / 

sm[ii N rep F (p)\ 


e -w( P ) 

Wp) 


l\ G o{p) + ^G 0 (p ) - i)G»{p) 


sin[7uF(p)] 


( 110 ) 


This formula is factorized into three essential parts. It contains a term e _H A p) jG"(p) which 
is responsible for an exponential decay of the probability amplitude of ionization v 4^ addle ' 1 (p) 
while increasing the photoelectron energy. Another factor, sin[|Go(p) + \Gq{p) — ipG"(p )\, 
corresponds to slow modulations of the probability amplitude on the electron energy scale. 
This is in contrast to the last term in Eq. (IllOj) . sin[7rlV re pF(j>)]/ sin[7rF(p)], which we call 
the diffraction term. As we are going to illustrate, this term is a source of very sharp 
peaks in the energy spectrum of photoelectrons, similar to diffraction fringes observed in the 
experiment by Davisson and Germer [H^ ]. 


21 



























10 12 14 16 18 



10 12 14 16 18 


Ep/ 

FIG. 6. (Color online) In the top panel, we show a portion of the energy spectrum presented in 
the upper frame of Fig. O Only the curve for N rep = 3 is plotted, and the results are not scaled by 
7V r 2 e p. Vertical lines mark the energies at which we observe the main maxima. The same but for 
IV r ep = 10 is plotted in the middle panel. Note that in both cases the main maxima occur at the 
exact same energies of the final electron. At those energies the function F ( p ), drawn in the bottom 
panel as the solid blue line, takes on integer values. Note that F(p ), in contrast to the dashed 
black line, is not a linear function of its argument. Therefore, the peaks in the energy distribution 
of photoelectrons are not equally spaced. 


2. Numerical illustrations 


In Fig. Owe plot the quantity |p||«4k-(p)| 2 which, according to Eq. (KJTjh is proportional 
to a triply differential probability distribution of ionization, 


\p\\A K (p)\ 2 ~ 


d 3 V K 

dQpdEp 


( 111 ) 


For a visual purpose, we have multiplied this distribution by . While this distribution 

is invariant with respect to the azimuthal angle, for the polar angle we have chosen 6 p = 0 . 27 T. 
In the upper frame, we present the exact results based on a direct numerical calculation of 
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the integral in Eq. (l96j) . The mirror-reflected curves, shown in the lower frame, have been 
calculated using the saddle point method with respect to the aforementioned integral, i.e., 
based on Eq. (I103jl . While spectra in both frames differ in magnitude, their actual patterns 
are the same. In each frame we present three curves. The solid black envelopes correspond 
to the case when the driving pulse is a single-cycle pulse (A rep = 1). As it follows from the 
saddle point treatment (II101) . in this case the diffraction term equals 1 and, therefore, only 
slow modulations of the spectra are manifested. As we have also checked, for more energetic 
photoelectrons we observe similar modulations which, however, decrease in magnitude. Such 
a behavior can be explained by the exponentially decaying term in Eq. (II 10|) . The dashed 
red line is for A rep = 2, meaning that the driving pulse consists of two one-cycle pulses. 
Already in this case, a diffraction pattern is observed. We see a very sharp peaks within the 
envelope. This happens for any N rep ^ 2, in agreement with formula (1 11 01) . For instance, for 
Arep = 3, the corresponding sharp peaks are plotted with the solid blue line. Note that each 
spectrum was divided by A r 2 ep . This resulted in nearly same heights of the peaks for different 
A rep . While for energies 6 cul < E p < 25cul, one can actually see that the scaled peaks have 
the same heights for different A rep , for energies 3 cul < E p < 6co^ this is not exactly the case. 
Such a behavior of the presented spectra can be explained using the derivation based on the 
saddle point approximation (1 110 [) . According to this formula, the sharp peaks appear at 
electron energies such that F(p) = L, where L is integer. This behavior is distorted by the 
term sin[|G 0 (p) + |GoCp) — ' l l J G"(p )\, which manifests strongly for 3 cul ;$ E p < 6cul- Let us 
also note that the spectra divided by A r 2 ep have contact points at such electron energies E p 
that the phase of the probability amplitude of ionization takes the same values regardless of 
N 

1 y rep • 

In the top panel of Fig. El we plot a portion of the spectrum presented in Fig. El for 
A rep = 3. The same but for N rep = 10 is plotted in the middle panel. Note that in both cases 
we observe the enhancement of the spectra at the exact same electron energies, as indicated 
by the solid vertical lines. In the bottom panel, we show the function F(p) (solid blue line). 
As expected, the main maxima in the upper panels occur at those photoelectron energies 
when F{p) takes integer values. At these energies, the diffraction term in (II101) tends to 
Arep and, hence, the respective probability distributions scale as A 2 cp . Since the major peaks 
become more narrow with increasing A rep , the angle-resolved probability of ionization, when 
integrated over the electron energy, scales approximately as A rep . Therefore, for the Keldysh 
theory it is meaningful to talk about the probability rate of ionization per one modulation 
of the laser pulse. Also, note that F(p) is not a linear function of the photoelectron kinetic 
energy, which is in contrast to a straight line (dashed black line) shown in the bottom panel 
as well. It means that, in general, the enhancement peaks are not equally spaced on the 
photoelectron energy scale. Another feature which can be observed in Figs. El and El is that 
with increasing A rep , there appear (A rep —2) additional maxima between any two consecutive 
main peaks. Their positions can be derived from Eq. (II101) . F{p) = L + (M + l/2)/N rep 
where M = 1, 2,..., A rep — 2. These additional maxima are accompanied by zeros in the 
energy spectra. For N vep ^ 2, there is always (A rep — 1) zeros which are observed when 
F(p) = L + M/N rep with M = 1,2,..., A rep - 1. 

The diffraction pattern in the photoelectron energy spectra is observed only when A rep ^ 
2, i.e., when the driving pulse is composed of at least two modulations. Its features can be 
explained based on an approximate formula for the probability amplitude of ionization (II 10)1 . 
which suggests a very intuitive interpretation of the observed pattern. Namely, the proba¬ 
bility amplitudes from each modulation interfere constructively, leading to enhancements at 
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certain electron energies. One can conclude, therefore, that each modulation acts as a slit in 


the Young-type experiment of matter waves performed by Davisson and Gerrner 57]. Note 


that similar diffraction patterns can be observed in other strong-field processes as well, with 
the most recent examples in the area of strong-held quantum and classical elcctrodynam- 
58 611 or in optics for electromagnetic waves passing through diffraction gratings j2|. 
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C. Combs in the GEA 


We have shown in the previous Section that diffraction patterns in the photoelectron en¬ 
ergy spectra follow from the Keldysh theory. Since the Keldysh theory neglects the Coulomb 
interaction between the ejected electron and the residual ion, the question arises whether the 
similar patterns can be still observed if the Coulomb interaction between the two is taken 
into account. To answer this question we will use now the GEA. 

The probability amplitude of ionization in the first order of eikonal perturbation the¬ 
ory 1)58]) . contains the extra space- and time-dependent phase factor as compared to the 
Keldysh amplitude (|9T]) . This factor functionally depends on the real-time classical trajec¬ 
tory in the laser field, i.e., 


y^o(p) r 2 ™ 


A {1 \p) = - i-- / d0e iG(p ’ < « / dVe- iqWV (-e^(</>) ■ r')Mr')e~ iU[r '^ plrcl \ (112) 


u 


where 


U[r’,(j),p\r c i] = 


and [cf. Eq. 
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r d (a; r',0,p) = R p (r',-, 
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Note that the real-time classical trajectory r c i(a; r', 0,p) depends on r' and p (here a plays 
the role of time in units of 1/cu) through the initial and final conditions, 


r cl (a]r',(j) : p) 


- 7 ^r c i(a] r, 0, p) 


= r i 

P 


ct=2tt mu 


(115) 

(116) 


respectively. It is worth noting that, by following the standard procedure (see, e.g., Ref. 
the functional 

W[(p, p\r cl } = -G(p , 2vr) + G(p , 0) - U[r', 0, p\r cl ] (117) 

can be rewritten in the form 


W[^ p|r d ] = S(i)>|r d ] + I»r d ((() • r’ a (<p) - p ■ r d (2ir), 


(118) 


where we have used the abbreviation r c i(cr) = r c\{o",r',(j),p). Here, the ‘prime’ means the 
derivative over a and 


1 /‘ 27r 

S[(j)\r d ] = - da£ cff (r cl (a),r' 1 (cr),cr), (119) 

w J<f> 
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is the classical action with the effective Lagrangian 


C e ft{r cl (a),r' cl ((j),(j) = ^-[r' d (cr )] 2 + eS(a) ■ r d (a) - (r A (a), a ~ 4>) + E 0 . (120) 

Here, the effective potential is defined by Eq. (148|) with time in units of 1/co. Since we deal 
with the static Coulomb potential then depends only on the classical trajectory and 
the phase difference a — 0. Finally, the probability amplitude (1 112 j) can be put in the form 

i(p 2 /2m-E 0 )T r2n r 

A {l \p) = -i- / d0 / dVe- iqWV (-e£(</>) ■ r')%(r')e lW ^ p ^\ (121) 

u Jo J 

which is suitable for the saddle point and quantum trajectory analysis. 


1. Quantum trajectories 

The direct integration over the space variables, d 3 r', in (1 121 j) is very difficult to carry 
out as the effective potential, being the function of a complex argument, oscillates rapidly. 
For this reason, it is very convenient to apply the saddle point method. In our further 
analysis we use the simplest approximation, namely, in the effective potential we replace 
the classical real-time trajectory r c Ja) by its quantum analog (which is frequently called 
the complex-time trajectory |38|) being the solution of the free particle Newton equation 
in a laser held. This is in agreement with the assumption, which is commonly made in 
the strong-held approximation, that the binding potential rather marginally modifies the 
electron trajectory in the laser held. 

The quantum trajectory is the solution of the classical Newton equation in the laser held, 

( 122 ) 

H mu z 

where the ‘ prime’ again means the derivative with respect to the phase a. These trajectories, 
however, have to fulfill the complex initial conditions for a = <p s , 

Re[r q (<j) s ;p, <j> a )} = 0 , [,mu}r' q ((j ) s ; p, (j) s )] 2 = -A 2 , (123) 

where A relates to the binding energy of a hydrogen-like atom (j92]h Following Ref. j38j, we 
can write down the quantum trajectories of the form 

r q (a]p,(j) s ) = — [cr - Re(0 s )] + — ct(a) - —Re[a(<&,)], (124) 

muj m m 

where [cf., Eq. f[79jl ] 

«(a) = —— [ A(ip)dip = -lf a {cr)e z . (125) 

OJ Jq uj 

These trajectories satisfy both conditions (II23j) . The hrst condition in (j!23j) states that the 
real part of the complex trajectory starts at the center of the atom. The second condition 
defines the initial phase (f> s , which turns out to be the saddle point solution of Eq. (|99]h In 
addition, the trajectories (1 124ft are real in real phase a, 

Im[r g (Re a; p, <j) s )} = 0. (126) 
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FIG. 7. (Color online) The same as in Fig. [5} except that in the lower frame the spectra have been 
calculated based on the GEA with the application of the quantum trajectory method [Eq. (11311) ]. 


The same holds for the velocity, 


Im[ r g(Re P, 4>s)} = 0. (127) 

Since we deal with the Coulomb-free trajectories, if the laser pulse is switched off the pho¬ 
toelectron will carry the momentum p. In other words, 

mur' q ( 2vr; p, <j> 3 ) = p. (128) 

This condition agrees with the assumption made above that the electron final state is ap¬ 
proximated by the plane wave. It also agrees with the numerical analysis showing that the 
quantum trajectories at the end of sufficiently intense laser pulses (which is the case for the 
intensity considered in this paper) are far away from the Coulomb center. In other words, 
that the electron final momentum is rather marginally affected by the interaction with the 
residual ion. 

For each saddle point (p s and the corresponding quantum trajectory, we can define the 
generalized eikonal 


Xq\pAs) = 
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and its original counterpart 


Xg ,original (P; 0s) 


Zac 

u 



da 


1 

r q (<r, P, 


(130) 


With these definitions the probability amplitude in the first order of eikonal perturbation 
theory and in the saddle-point approximation adopts the form [cf., Eq. (11031) 


•^saddle (p) 



^rA \ 2 

mu) 6 ^ G"(p,(j ) s ) 


(131) 
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FIG. 8. (Color online) The thin black line represents the energy distribution of ionization in the 
saddle point Keldysh approximation, Eq. (1103)1 . and the thick blue line corresponds to the GEA, 
Eq. (1131 j) . The dashed red line depicts the result for the EA. The upper and the lower panels show 
the distributions for A r rep = 1 or 3, respectively. The remaining laser pulse parameters are the 
same as in Fig. [5j 


and similarly for the EA, with the replacement of Xq^ by Xq, original- 

In Fig. [71 we present the energy distributions of ionized electrons similar to Fig. El The 
difference is that, this time, the lower frame shows the spectra calculated within the GEA 
and saddle point approximation (1 131 j) . In this case, we account for the Coulomb interaction 
between ejected photoelectrons and their parent ions. As we see in the lower frame, the 
positions of peaks and zeros in the spectra are almost identical as in the upper frame where we 
plot the spectra calculated based on the Keldysh approach (1961) . What is changed, however, 
when we account for the Coulomb interaction between the electrons and the residual ions, 
is a significant enhancement of the ionization signal. Also, we observe a partial loss of 
coherence since the distributions do not scale any longer like iV r 2 ep . 

The enhancement of the ionization yield can create some doubts about the validity of 
eikonal perturbation theory. Let us note, however, that the perturbation is carried out in the 
exponent. In this particular case, the applicability condition for this approximation is such 
that in Eq. (I13ip the eikonal term X^ijPi^s) should be much smaller than the zeroth-order 
term, G(p,(j) s ), for both the real and imaginary parts. This condition is very well fulfilled 
for the laser pulse intensity considered in this paper. 
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2. Prospects for using the GEA 


Now the question arises: To what extent the GEA is better than the EA? To answer 
this question, in Fig. [8] we compare the predictions of both approaches. The thin black 
line represents the results calculated based on the Keldysh theory while the blue thick line 
is for the GEA. The EA results are represented by the red dashed line. In each case, the 
saddle point method was used. The energy spectra presented in the upper panel are for 
the one-cycle pulse, N Tev = 1. We observe a rather marginal difference between the GEA 
and EA results, and a significant enhancement (by roughly one order of magnitude) of 
these distributions as compared to the Keldysh approximation. Qualitatively, however, all 
three distributions look similar. The differences appear for longer pulses, when N iep > 1. 
This is illustrated in the lower panel for N rep = 3. For electron kinetic energy E p fa 3o;l the 
distributions for generalized and original eikonals differ significantly. However, for E p > 5o;l 
(which is not shown in the figure) both approaches again give nearly the same results. The 
wiggles observed for the EA distribution can be explained if we note that for the laser 
field parameters considered in this figure the ponderomotive energy is close to which 
means that the structure appears for E p fa 3 U p . It is well-known that for such energetic 
photoelectrons some of the complex trajectories can return very close to the origin of the 
Coulomb potential. As we have checked, this is the case here. Since the original eikonal is 
singular for such trajectories, we observe the rapid change of Xq, original (p, 4>s) when the kinetic 
energy passes through the value 3U P (in the considered case the real parts of x q , original (p, 4> s ) 
exhibit the sharp peaks for these particular trajectories). This results in wiggles observed 
in the lower panel of Fig. [S] for the EA. Such a behavior, however, is not observed for the 
GEA, as it is not singular for trajectories returning to the potential origin. If we compare 
the GEA with the Keldysh approach, we see the enhancement of ionization but again the 
distributions are qualitatively similar. In our opinion, the lack of spurious behavior for 
trajectories returning back to the vicinity of the parent ion and the fact that the first Born 
approximation is the limiting case of the GEA make the approach presented in this paper 
an attractive tool for investigations of ionization, rescattering, and high-order harmonic 
generation by strong laser pulses. This includes also more complex systems such as two- 
atom molecules or fullcrenes. 

In the literature (see, e.g., Ref. [33j) two names for the method, quantum trajectories 
and complex-time trajectories, are frequently used. In light of our analysis and more thor¬ 
ough studies carried out, for instance, in Refs. 33], 38], the second name seems to be more 


appropriate, as the trajectory r q {a\p,(f) s ) satisfies the Newton equation with the classical 
binding potential V(r) but with the complex initial conditions. In such a formulation of the 
method, there are no quantum signatures in the definition of r q {cr,p,(j) s ). This approach, 
however, leads to some problems related to the Coulomb singularity at the origin. In quan¬ 
tum theory (for the Schrodinger and Dirac equations), this singularity does not create any 
difficulties. Our investigations show that this obstacle can also be eliminated in the complex¬ 
time method. Indeed, the form of the effective Lagrangian (I120p suggests to assume that, 


up to the first order of eikonal perturbation theory, the trajectories should fulfill the Newton 
equation of the form [cf., Eq. I1122j) ] . 


muj 2 r"(a-,p,(j ) s ) = e£(a) - Wi 1 ) (r,(ff;p,^),ff - «#.»), 


(132) 


with a suitable initial conditions. The effective potential !()) '( r. o). contrary to the classical 
one V(r), is not singular at the origin for non-zero time and is smeared out by the ‘ quantum 











diffusion ’ represented by the Laplacian and the non-linear term in Eqs. (126|) or (1281) . In 
other words, it accounts for the spreading of the electron wave packet during the quantum 
time evolution. The Laplacian introduces the Planck constant into the definition of the 
effective potential in the first order eikonal perturbation theory. In other words, \r, a) 
differs from V{r ) by quantum corrections which vanish in the limit h —> 0. This also means 
that r q (a] p, (f s ) does contain quantum corrections and, therefore, we should rather call 
those trajectories ‘ complex-time quantum trajectories , . The effects related to the quantum 
corrections in r q (a-,p,(j) s ) are now under investigations. 


IV. CONCLUSIONS 

We have formulated the GEA for ionization processes driven by strong laser pulses. As 
we have shown, the Born approximation arises as the limiting case of our approach. The EA 
does not have this property, which significantly diminishes its applicability to the rescattering 
phenomena. Moreover, the EA is singular for the trajectories that come back to the center 
of atomic potential. We have demonstrated that the GEA does not have this shortcoming 
either. This makes it a very promising tool to study rescattering-related phenomena, with 
the most prominent example of high-order harmonic generation. 

Using the GEA we have discussed the appearance of coherent diffraction patterns in 
photoelectron energy spectra and their modifications induced by the interaction of photo¬ 
electrons with the parent ion. We have identified the conditions necessary to obtain such 
coherent patterns. If a pulse consists of at least two modulations, each of these modula¬ 
tions acts as a slit in the Young-type experiment for matter waves resulting in a coherent 
enhancement of ionization signal at particular electron energies. As we have illustrated this 
numerically, if we increase the number of modulations within a pulse, the comb-like struc¬ 
tures in the energy spectrum of photoelectrons become similar to the Alike structures. This 
is particularly interesting in the context of designing new sources of electron pulses, which 
is another topic to be studied in near future. 
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